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Abstract 

Quasi-two-dimensional transport is investigated in a system consisting of one ferromagnetic layer 
placed between two insulating layers. Using the mechanism of skew-scattering to describe the 
Extraordinary Hall Effect (EHE) and calculating the conductivity tensor, we compare the quasi- 
two-dimensional Hall resistance with the resistance of a massive sample. In this study a new 
mechanism of EHE (geometric mechanism of EHE) due to non-ideal interfaces and volume defects 
is also proposed. 

PACS numbers: 75.47.-m, 75.70.-i, 73.50.-b 
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I. INTRODUCTION 



Recently there has been an increased focus on the fabrication of new spintronic devices 
based on the tunnel magnetoresistance effect (TMR)^ and the recently discovered spin torque 
effect (ST)2 in multilayered structures consisting of several ferromagnetic nanolayers sepa- 
rated by thin insulating barriers. The geometric parameters of such structures also give us 
the possibility to investigate some of their quasi-two-dimensional transport properties and 
particularly to compare the relationship between diagonal and off-diagonal (responsible for 
EHE) conductivities in massive and above mentioned samples. This study is of great inter- 
est to the field of spintronics since quasi-two-dimensional EHE may provide an additional 
mechanism for recording and storing information in logic devices. Focusing on this aspect of 
EHE we study both diagonal and off-diagonal (Hall) conductivities in a system consisting of 
one ferromagnetic layer of thickness a placed between two insulating layers and magnetized 
in the direction perpendicular to the interfaces (z— direction). 

The remainder of the paper is organized as follows. In Section II EHE is attributed to the 
mechanism of skew-scattering^ on the bulk and interface impurities. We report our results 
for the conductivity tensor including size effect terms calculated within the framework of the 
Kubo formalism and compare the quasi-two-dimensional Hall resistance p^o with that of a 
massive sample p^ ;fc . Section III discusses the nonideality of the interfaces and the existence 
of volume defects that influence the form of current lines. We propose a new mechanism 
of EHE (we will refer to it as the geometric mechanism of EHE) due to these defects using 
the diffusion equation^ and taking into account that the diffusion coefficient has off-diagonal 
component proportional to the spin-orbit interaction. We summarize our results and offer 
some conclusions in Section IV. 



II. SKEW-SCATTERING MECHANISM OF EHE IN A THREE-LAYERED 
STRUCTURE 

We will consider the geometry of current parallel to the y— direction resulting in the 
appearance of the x— component of the Hall field. In this case: 

3y = a yyE y + &yxE x j x = a xy E y + a xx E^ = (1) 

t?h _ ° x y jp „h _ a yx (c)\ 

x — n- y H2D,bulk ~ , 2 V > 

where a a p are diagonal and off-diagonal conductivities, E^ is the Hall field. For the calcu- 
lation of p2 D bu i k we will use the Kubo formula with vertex corrections responsible for the 
transverse component of the current: 

he 2 \ - 

= —4 22 u ^'p g kA zz ') g k> k ( z ' z ) ( 3 ) 

nap. * — ' 

kk'z' 
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where v K is the velocity vector along the interface, Gq k (zz'), Gq k (zz') are advanced and 
retarded Green's functions in mixed coordinate-momentum representation. To calculate the 
y — x component of the conductivity tensor we will use the perturbation theory and will 
take into account only the corrections in linear order of the spin-orbit interaction: 

GlAzz') = G+ K (zz')5 KK , + G+ (zz") J2 ( T kAp z ") + KAp z ")) (4) 

p 

H S K °APZ") = (PZ") 4™* [« >< «'L = W) <%™>z («*«!, - <Ky) (5) 

where T^ K ,{pz) and H^° K ,(pz) are the scattering matrix and the spin-orbit interaction, cor- 
respondently, dependant on the type of atom in (pz) position ; m z is the unit vector along 
the magnetization, \™ is the spin-orbit parameter and ao is the lattice constant. It follows 
from Eq. © that H%(pz") = -H$ K (pz"). 

For the T— matrix in one-site approximation we can write down: 

T KK ,(PZ) = E e^'^ \ %7lr\- - ^ ® 

^ 1- (e n - E (z)) G (p n z, p n z) 

K° K ' (P z ) = E e^'^-^\Z (7) 

rn 

where e n is the one-site energy. For the binary system AB e n and A^ take values €a,b and 

2 -^2lL 

X a A ° B , respectively; G (p^z, p^z) = ^J °° kcIkG k (zz). £ is the coherent potential and can 
be found from the system of self-consistent equations. The first equation of the system is 
valid for scattering on both bulk and interface impurities: 

Ca £ A - g (g) + ^ e B - S (z) = Q 

1 - (e A - S (z)) G (p" n z,p~* n z) B l - (e b - £ (z)) Go (p" n z, p" n z) 

while the second one is written in the form corresponding to the interface scattring since we 
are interested in calculation of the interface coherent potential: 

n ( - -> \ G (p n z,p n z) 

G {p n Z, pnZ) = ~ —, . _ — 9 

l-E(z) G (p n z,p n z) 

For the calculation of a xx , a yy we will use Eq. ([3D with Green's functions diagonal on k 
and renormalized on the coherent potential. For off-diagonal component averaging on the 
impurities distribution gives: 

O yx = E K ( ZZ ")\ 2 K (Z'z't 1 ™ ( T :AP^")K' K W)) (10) 

where p" n z" is the impurity position. We keep in Eq. ( TTOl ) only the main term with n = m. 

For the binary AB structure and purely random distribution of A, B summing over p^ n 
will give S KK ". It is convenient to divide A so into average and scattering parts: 

Aa = c a Xa + c b Xb + c B X A ° ~ c B \% = A + c B 5X so (11) 
A^° = c A X s B ° + c B X%° + c A X s A ° - c A X s A ° = A - c A SX so (12) 
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The average in Eq. (1101) is: 



^/ £ a - g (g) + ^ e_B - E (g) | ^ 

A l - (ca - S(2))G(pn2,p^) B l - (e B - E (z)) G{p n z,p n z) 



+ <5A so c A c B 



e A - E (s) e B - £ (z) 



1 - (e A - £ (z)) G (p~* n z, p* n z) 1 - (e B - E (2)) G (p~;z, pn^) 

According to Eq. and Eq. ([9]) we only take the imaginary part of the last term in Eq. 
13l) . We rewrite Eq. (fl~3l) with renormalized Green's function: 

im-|- £ (T B (z)ATW> - *A«Wm ^ ^ " S W) (1 " S W G ° ( ^' 



N „ n V 1 - e A G (p^z, p^) 

- £ (z)) (1 - £ (z) G (p~;z, pn^)) 



(14) 



1 - e B G (p" n z,pZz) 
G {p" n z, p" n z) = F (z) (15) 
Eq. (J6j) for coherent potential is: 

^ c A c B 5 2 F(E,z) 
E(£ "~ ) = £+ 1-(I-£( £ , 2 ))F( £ ,;) <16) 

e = c A e A + c B e B , I = c A e B + c B e Al 5 = e A -e B (17) 

Usually it is convenient to choose e = 0. In this case e = — (c A — c B ) 5, e A = c B 8, 
t B = -ca6. 

Next we assume that scattering parameters for bulk impurities such as scattering potential 
and concentration are small enough to keep only the main terms for all values. In this case 
the imaginary part of £ is of order S 2 so for a^ k we will use: 

(Tn{ z )^?( z )) ~ ^Kulk^bulk c AbulkC B bulk {c A bulk ~ CBbulk) 

! (18) 

x ImF 2 2 

(1 — c Bbu i k S bu i k ReF) (1 + c Abu i k S bu i k ReF) 

For the interface values the full self-consistent scheme is necessary. For both bulk and 
interface we suppose that the real part of the coherent potential just represents the renor- 
malization of electron spectrum so E can be considered as purely imaginary. 

The zero order Green's function found from Schrodinger equation in k — z representation 
is (we will further use the units with energy dimension [L\ = A): 

1 



Gq k (0 < z' < z < zi) 



2«jfci (e ikia (q + ikif - e~ ik ^ a (q - ih) 2 ) 
x ( e ikl{z - Zl) (q - ih) - e - ikl{z ' Zl) (q + ih)) ( 19 ) 
x (e lkiz ' {q + ih) - e~ iklZ ' (q - ifc x )) 
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where and z\ are the coordinates of the left and right interfaces, a is the layer thickness; k F , 
li are the Fermi momentum and the mean free path for spin "up", respectively, C\d\ = -f 1 
(for spin "down" we will use index 2); U is the height of the potential barrier. 

The poles of the Green's function in Eq. (fT9l) define the quantized energy spectrum of 
the thin ferromagnetic layer. 



A. Calculation of the bulk quasi-two-dimensional diagonal conductivity 

For a xx = o yy in Eq. ([3j) we will take into account scattering on the interface respon- 
sible for size effect as well as on the bulk of the sample by using the Dyson equation with 
renormalized Green's function: 

G K (zz') = G 0K (zz') + G 0K (zO) SG k (0/) = GW (zz') + ^(^)SG 0k (0^) (24) 

Integrating over z' from to z for z' < z and from z to a for z' > z gives the conductivity 
in the units ohm~ x cmT x : 

t o-o^ilO 8 f K z dnNom} . . 

aL = ^rJ^D^- <25) 

jVom T = (q 2 + cj) [(q 2 + c{ + |£| 2 + 2qReT) sinh 2^ + 2ci \ImE\ cosh2d 1 a] 

+ 2c\lmYj~ [(g 2 + c 2 ) cosh 2c?! (z — a) + 2 sinh 2 ,2 (cos2cx (z — a) + 2gcisin2ci (z — a))] 
- 2ci (q + ReT.) (q 2 + c 2 ) sinh 2d x {z - a) sin 2c x z + (q 2 + c\ + |£| 2 + 2g J Re£) 
x [(g 2 + c 2 ) sinh2di (z — a) cos2ci2; — siah2diz ((g 2 — c 2 ) cos2ci (z — a) + 2gcisin2ci (z — a))] 

(26) 

Den T = (q 2 + c\) (((q 2 + cf) + |S| 2 + 2g J ReS) cosh 2d x a + 2c x |/mE| sinh 2d x a) 

- [{q 4 ~ Qq 2 cl + 4) + (|S| 2 + 2qReY J ) (q 2 - c 2 ) - Aqc{ReY] cos2 Cl a (27) 

-2qd (2 (q 2 -c\ + qReT) + |S| 2 ) sin2cia 
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Figure 1: Averaged diagonal conductivity as a function of a (thickness): < o xx > (a) for k F 
l.liA- 1 ), kp = (^(i- 1 ), h = 100(A), h = 60(A), c = 0.3 (see Eq. ([Mil ) 



where <t = ^ = ^(o/im -1 ) is the elementary conductivity of one channel. 

For large enough layer thickness we can average Eq. ( !25|) over oscillations so that the 
averaged conductivity is: 



2ixk\ 



Cl 



h 




Ci sinh 2dxa 


a (q 2 + cf + 


E 


2 + 2gi2eS) sinh2dia + 2c x 


ImTi 


cosh2<iia 



(28) 

The first term in Eq. ([281 is the conductivity of the massive sample and the second one 
is due to the quasi-classical size effect. The full conductivity representing the sum of two 
spin channels is shown on Fig. (pQ) as a function of the layer thickness. 



B. Calculation of the off-diagonal conductivity due to the spin-orbit interface scat- 
tering 



Now we will calculate &l y (z) using Eq. ( fTOl) with Green's function defined by Eq. ([24 
and z" = 0. After integration over z' the conductivity is: 



a 



■HI 



2>K 2 k\ 



k 3 (Ikci 



Im (T (p n Q) X so (p n 0)) (q 2 + c\) sinh 2d 1 a 



D ikia 



x / k ■ dn 



3 7 (<? 2 + c i) cosri 2c?i (-2 — a) — (q 2 — c 2 ) cos 2ci (2 — a) — 2cig sin 2ci (2 — a) 



(29) 
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Figure 2: Off-diagonal conductivity as a function of a (thickness): a xy (a) for k F = 1.1 (A 1 ), 
kp = 0.6(i- 1 ), h = 100(A), Z 2 = 60(A), c = 0.3, A so = 0.05(i4 _1 ) (see Eq. (pjl ) 

This conductivity oscillates with the thickness and the distance from the interface z = 0. 
Its behavior becomes more clear after averaging over oscillations: 

/ I v = ap/i^lO 8 /" Im (T (p»0) A*° (p n 0)) sinh 2rf lQ 

\ a W S7r2k y J K KC i ^ 2 + c 2 + | S _|2 + 2q ReJ:-) sinh 2d 1 a + 2c 1 /TO£- cosh 2d x a 

/, , cosh 2d\ (z — a) 
K UjKj 
(g 2 + c\ + |S-| 2 + 2qReL-) sinh2rfia + 2c 1 JmE" cosh2dia 

(30) 

The same is done for spin "down". The sum of these two terms is shown on Fig. ((2j). 

It is clear that a\ y decreases with z — ► a since the functions cosh 2d±(z— a) have maximum 
values at z = 0. Averaging o\j^ over these functions gives the factor ^ so for infinite a this 
term tends to zero. 

We also calculate the bulk conductivity a\ y + a^. y (see Fig. ((2D) using Eq. (TTOl with 
additional integration over z" and bulk scattering parameters with the bulk coherent poten- 
tial in Bohrn approximation T, bulk = icbuikiA ~ c buik)^ u i k ImF bulk (zz). In the absence of the 
interfacial scattering this approach gives us: 
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Figure 3: Bulk off-diagonal conductivity as a function of a (thickness): a b y lk (a) for k F = LI (A 1 ), 
k l F = O^A" 1 ), h = 100(A), h = 60(A), X s b ° ulk = 0.03(A), c bulk = 0.01, S Mk = l^A- 1 ) (see Eq. 



xy 



a l 2 a 3 Im {T bulk \ s b ° ulk ) 10 s y (g 2 + c\f sinh 2^a ^ f t^dn 1 2 2 . 2 



TJ*J* = __ . , _^ ~ _ x 



87r 2 ^ 2 7 ci Deri J c x Dev} 



x sinh 2d\a — (g 4 — c 4 ) (sinh 2rfi,s cos 2ci (z — a) — sinh 2di (z — a) cos 2c\z) 
+2qc\ (sinh2di2;sin2cx {z — a) + sinh2<ii (z — a) sin2ci2;)} 

(31) 

.Den 1 ^ = cosh2c/ia (g 2 + c 2 ) 2 — (g 4 — 6g 2 c 2 + c 4 ) cos2cia + Aciq (g 2 — c 2 ) sin2cia (32) 
or after averaging over oscillations: 

This conductivity has an oscillating behavior for the thin layer but tends to the constant 
value when a — > oo which coincides with its value for the massive sample. If we take into 
account the interfacial scattering the expression for a b y lk becomes too complicated so we 
don't show it here. But the thickness dependences of a^ k and — = — calculated using 

r X V p Vxx ° 

the full formula with renormalized Green's function are presented at Fig. (plj and Fig. j5)), 
correspondently. The bulk parameters are: X s b ° ulk = 0.03(A- 1 ), c bulk = 0.01, 5 bulk = 1(A" 1 ). 
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Figure 4: Bulk off-diagonal conductivity with interfacial scattering as a function of a (thickness): 
a h ^ k {A) for k\ = l-UA- 1 ), k l F = O^i" 1 ), h = 100(i), l 2 = 60(A) 
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Figure 5: Hall angle ^- = as a function of a (thickness); bulk parameters: \'l°i k = 0.03( J 4 _1 ), 
c buik = o oi ; s hulk = 1{A~ 1 )- interface parameters: A so = O.OS^- 1 ), c = 0.3, 5 = 1.5(A _1 ) 

III. GEOMETRIC MECHANISM OF EHE 

Let us consider an electric current through a thin ferromagnetic layer of thickness a 
located between two thin insulating barriers and magnetized in z— direction perpendicular 
to the interfaces. The interfaces are not ideal and besides the impurities have topological 
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defects which will be modeled as cylinders of radius R so that the current lines in the vicinity 
of these defects follow their shape. 

Diffusion equations for charge and spin currents in the absence of precession are: 

dn dj x djy 

dt ~*~ dx dy ^ 

?™L + ?2k + ®2k = -™L (35) 
dt dx dy r a f 

For the stable state solution || = ^ = 0. For currents we have the system of equations: 

^ ^ dn ^ „,dm _ dn ^ ^,5m , , 

f e = o xx E - D xx — - D xx (3>— - D xy — - D xy p>— (36) 

dfi dm dft dm, 

fe = °yy^ + o^E - D yy — - D yy (3'— - D yx — - D yx (3'— (37) 

f m = Po xx E - D xx p'- - D xx — - D xy $- - D xy — (38) 

f/71 (fill f/71 (/Til 

Jl = P*yy0 + ^ yx E - D w ff— - D yy — - D yx /3'— - D yx — (39) 

Here we take into account that E = {E 1 , 0,0}, m = {0,0, m} are the electric field and 
spin accumulation correspondently, D a/3 are the components of diffusion coefficient tensors. 
Off-diagonal components D xy and D yx of these tensors are proportional to the spin-orbit 
interaction and they are antisymmetrical in the x — y transposition. For a metal with the 
cubic symmetry o xx = a yy , D xx = D yy = D . Then we insert Eq. (1361) -(1391) into Eq. ( |34l ) 
and Eq. ( 1351) and after some manipulations we obtain two equations: 

An = -(3'Am (40) 

771 

P'An + Am = — (41) 



And for T sf D (l - (5 12 ) = A 



2 . 



rn 



Am-- r = (42) 
X sf 

For the cylindrical defect shape it is convenient to search for solution in polar coordinates 
so we can rewrite Eq. (jl2l) : 

1 d ( dm\ 1 d 2 m m _ d 2 m 1 dm 1 d 2 m m 
r dr \ dr J r 2 dtp 2 \ 2 j dr 2 r dr r 2 dip 2 X 2 s j 

where ip is the angle between x— axe and the radius- vector r with the coordinates (x,y). 
The solution of Eq. @HJ) is: 

m = mi (r) m 2 ((f) (44) 
m 2 (</?) = A\ n cos n<p + A 2 n sin nip (45) 
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As = —m 2 n 2 , Eq. ( l43l) can be transformed: 



m 2 (y?) 



d 2 rrii 1 <9toi 



n 

-z I mi 



The solution of Eq. (1461) is 5 : 



mi (r) = SfcKfc ( — 



m 



(^in cos n( P + ^2«. sin ny>) K k ( — 

\A s f 



where Kk(-c-) is the solution of the modified Bessel equation 5 . 



From Eq. ([40]) it follows that: 



n = —0m + uq 
An = 



For no the solution is: 



n 



^ (C ln cos ny? + C 2 „ sin ny?) — 



(46) 



(47) 
(48) 



(49) 
(50) 



(51) 



Taking into account Eq. ( l40l) we can rewrite Eq. (l36l ) and Eq. (l37fl : 

■a: _ rp n ^0 „ 9n 
Je — O^-fr - U 



j y — (T E — D 



tin 



' dx 
drip 
dy 



+ D 



: y dy 
dx 



,y.'- 



(52) 
(53) 



Now it is convenient to use the polar coordinate system and to write down r— and 
99—projections of the currents. Then we can use the boundary conditions to find unknown 
coefficients. These projection are: 



jre = °xxE cos (p + o xy E sin if 



(usual term) 



r. n dn dn dn . dn . 

dj re = ~D XX — COS tp - D xy~Q^ COS f ~ D yV~£fy~ Sm & + XV ~d~X ^ ^ 

(additional diffusion term) 

Now we will make some transformations: 

dn dn dr dn &Lp 

dx dr dx dip dx 

dn dn dr dn dip 

dy dr dy dip dy 



(54) 



(55) 



(56) 
(57) 
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Using the expressions for derivatives of r, ip over x, y, which is not too difficult to obtain, 



we write down the charge and spin currents in polar coordinates: 

Sj re = -D 



-D 



dn dr dn dip \ . / dn dr dn dip 

dr dx dip dx ) ^ \dr dy dip dy 

{ dn dr dn dip \ ( dn dr dn dip 

^ \dr dx dip dx) ^ \dr dy dip dy 



dn dn 

— D()— D X y—— 

dr rdip 



(58) 



5j 



m 
rrn 



D ?f - D^J- *L 
dr rdtp 

dm dm 

Do— D xy —— 

dr rdip 



Sjrm = P(T XX E COS If + (3(T yx E Sm ip 



(59) 

(60) 
(61) 



To find the unknown coefficients in Eq. (1481) and Eq. (|51~T) we will use the boundary 
conditions on the surface of the cylinder representing that r— projection of currents are 
equal to zero: 



r RO +sr R = o^D ^ R +D x , (M < 



n 



dr ■ - *V R Q \r=R = 'Do J2 ( Cl " COS ™f> + C m Sil1 U ^ R n+1 

n 

Dxy ^ (Cin sin nip - C 2n cos rup) 



R n+1 



It gives us the system 



with solution 



n 



D a. 



PxxE — —R 2 (DqC\i — D xy Ci\) 
<J xy E = R 2 (D xy Cn + D0C21) 



• ■ V D X y(J X y ^ Do<J X y ~\~ D X yO X 



D 2 o + Dly — D 2 + Dl y 

R 2 E (Dqo xx - D xy a xy ) cos ip - (Do<7 xy + D xy (Jxx) sin ip 
r D 2 + Dl 



(62) 



(63) 
(64) 



(65) 
(66) 



1 " xy 

Spin current r— projection is also zero, and we can write down: 

QtTI 0711/ 

j°R m + fijR = P E (a xx cos tp - a xy sin <p) = D — | r=iJ + D ^y R ^\r=R 

d ( r \ I f R 

= D (A u cos ip + A 2 i sin ip) — K x y~ J \r=R + D xy (-A n sin p + A 21 cos p) —K x y— 



(67) 
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Then we use some properties of Bessel functions^: 

K x (x) = lim — ] {I v _ x - I v+l ) = - lim -^—-I v (x) = (68 ) 

l siri7T [y + 1) ^^0siri7r^ x x 

8Ki (x) 7T 8 , , i , 7T 1v d , , . Ji (x) . . 

fl U = lim — t — TTT J *>-i - 4+1 = - I™ ttL (x) = 69 

ox ^-+Q I sin 7r [y + 1) ax k-+o sin itv x ox x 

here £ = y~, = t~- 

Kf ' or \ s f 

Inserting Eq. flEBD-flESD in Ec l- dSZj) we get: 



(3E (a xx cos tp - a xy sin tp) = -D (A u cos tp + A 2 i sin y?) ^ 

+ (A n sin <p - A 2 i cos tp) ^f- 



and it follows that: 



(3ER 2 [a xx D h i? + a xy D xy I { 



^ = — , v ; ":v , ' t : L 7 (73) 



A 



Dohi^j R +[D xy I ^)X sf 



m = — 



(3ER X s fl 



X <j rosj ( n, , /;•;,/, ( ^-J R + a xy D xy I ) A, 

+ sin I a xx D xy I I — ) A s/ - a xy D h ( — 1 R 



(70) 



= ~AAi| - D xy Aj-^f (71) 
/^o^ = -D^An^ - A)^2i | (72) 



f3ER 2 (a xx D xy I (^) A s/ - a xy D h (^-) /,') 

= ; ; — \ — 72 — ; ; — ; ( 74 ) 



(75) 



A// \ A sf, 

Now we can define the additional Hall field due to this cylindrical interface defect con- 
sidering that Hall electrodes are the surfaces with coordinates y = a and y = —a. This field 
is proportional to n(a) — n(—a), n = —f3'm + uq, r = .."i . After integrating over tp from 
to 7r for the left surface and from 7r to 2n for the right one we will have: 

/ \ ER 2 r sin 2 tp (D xy a xx + D a xy ) nER 2 (D xy a xx + D a xy ) 
Ma) - M -a) = 2—l dv a(Dl + D%) ™ 
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The second term due to —(3'm is: 



2(3(3' ER 2 r , sin2 <P i^D X yh fe) A s/ - a^ZVi (^) i?) 
w m \a) — n m (—a) = / dip 1 — - (77) 

At last, we have to multiply Eq. ( 1751 ) and Eq . ( 1771) by the concentrations of defects and 
electron charge. 

IV. CONCLUSION 

It was shown that due to the additional scattering of electrons on the defects of the 
metal-insulator interfaces the total conductance decreases. From Eq. ( 1281) it follows that 
for small values of the ratio j the bulk conductivity is completely suppressed and effective 
conductivity is proportional to the effective scattering length on the interfaces instead of the 
bulk mean free path. Hall conductivity, if we don't take into account the additional scew- 
scattering on the interface, decreases with decreasing the thickness of the ferromagnetic 
metallic layer. However the contribution to the Hall conductivity due to the additional 
scew-scattering on the interface increases. So the important characteristic of the considered 
device, Hall angle ^ = is larger for the thin ferromagnetic layer compared to the bulk 
layer. Besides that, the influence of insulator columns penetrating into the metallic layer 
may further increase the value of the Hall effect. 
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